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Abstract 

We give a necessary and sufficient condition for a homogeneous Markov process taking 
values in R n to enjoy the time-inversion property of degree a. The condition sets the shape 
for the semigroup densities of the process and allows to further extend the class of known 
processes satisfying the time-inversion property. As an application we recover the result of 
Watanabe in [27] for continuous and conservative Markov processes on R+. As new examples 
we generalize Dunkl processes and construct a matrix-valued process with jumps related to the 
Wishart process by a skew-product representation. 
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1 Introduction 

Let {(X t ,t > 0); (P x )xm n } be a homogeneous Markov process with semigroup densities (assumed 
to exist and to be twice differentiable in the space and time variables): 

P t (x,dy) =p t (x,y)dy. (1.1) 

For all x £ R" and some a > 0, the process {(t a Xi,t > 0); P^} is Markov and in general 
inhomogeneous. 

Definition 1.1. The process {(At, t > 0); P^} is said to enjoy the time-inversion property of degree 
a if the Markov process {(t a Xi, t > 0); P x } is homogeneous. 

The problems associated to time- inversion of Markov processes are closely related to the so-called 
dual processes in probabilistic potential theory (see [10] and [20]). A pair of dual processes is a pair 
of Markov processes whose resolvents are conjugate relative to some measure. The trajectories of 
these processes are connected by means of time-reversal. However, time-reversal in general violates 
homogeneity and changes the path properties of the process. In this paper, we are mainly concerned 
with the first issue, and in particular, whether one can characterize the semigroup densities of a 
Markov process such that homogeneity is preserved under time-inversion. 

Celebrated examples of Markov processes, known to enjoy the time-inversion property for a = 1, 
are Brownian motions with drift in M" and Bessel processes with drift (see [21] and [22]). Gallardo 
and Yor [13] recently worked out a sufficient condition on the semigroup densities for a Markov 
process to enjoy the time-inversion property. Their argument extended the class of processes to 
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processes with jumps such as the Dunkl process [26] and matrix- valued processes such as the Wishart 
process [3]. The aim of the paper is to find a necessary and sufficient condition and to provide some 
new examples. 

Section [2] contains the main theorem of the paper, which is proved in section [3] using straightfor- 
ward analytical arguments. Section |4] considers an application of the theorem to Markov processes 
on K.+ . The result is shown to be strong enough to entirely characterize the class of diffusion pro- 
cesses on R + that enjoy the time-inversion property and thus provide a different proof than that 
of Watanabe in [27j- Section [5] gives new examples of processes that enjoy the time-inversion prop- 
erty. We review how the generalized Dunkl process fits the requirements of the theorem and then 
introduce a matrix-valued process with jumps. The relation of the latter process to the Wishart 
process mimics the relation between the one-dimensional Dunkl and Bessel processes. 



2 Markov processes which enjoy time-inversion 

Fix x € K™. Recall that under P x , the process (t a Xi,t > 0) is Markov, inhomogeneous, with 
transitional probability densities qf}{z 1 y), (s < i; z,y S R n ), which satisfy the following relation: 



E, r 

where 



f(t a Xi) | s a X,_ =z = I dy f(y) £}{z,y) (2.1) 



qi* (a,b) = r™ ^ " 7 M • ( 2 - 2 ) 

Pi(x,-) 

The process (t a Xi , t > 0) is not uniquely defined from the knowledge of the semigroup densities 
p t (x,y). Actually, there exists at least one transformation that leaves the semigroup densities 
q^t(a,b) unchanged: Doob's /i-transform. 

Definition 2.1. Doob's /i-transform is the transformation 

T h :P x \^^^-e-^P x \ rt . (2.3) 
h(x) 

for some function h and some constant v > 0. 

This remark leads to our first assertion. 

Proposition 2.2. Two processes related by h-transforms yield the same process by time-inversion. 

It is hence legitimate to research for a criterium to classify all processes that enjoy the time- 
inversion property up to ^-transforms. The following theorem gives a concise statement of our main 
result: 

Theorem 2.3. The Markov process (t a Xi,t > 0) is homogeneous if and only if the semigroup 
densities of (X t , t > 0), assumed twice differentiable, are of the form: 

or if they are in h-transform relationship with it. The functions 0, p have the following properties 
for A > 0: 

1. $(\x,y) = $(x,\y); 

2. p(Xx) = A" p(x); 
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3. 9(Xy) = X 9{y) for some (3 € K; 

Moreover, if the symmetry condition = <!>(?/, x) is satisfied, then the semigroup densities are 

related as follows: 

g* (a0 (M) = J|^4 exp{tp(x)} Pt (a,b). (2.5) 

Such a decomposition of the semigroup densities shows furthermore that /i-transforms are the 
only transformations that leave semigroup densities of the time-inverted process (t a Xi , t > 0) 
unchanged. 

Definition 2.4. A Markov process (X t ,t > 0) is called semi-stable with index 7 in the sense of 
Lamperti [22j if: 

{(X ct ,t >0);P X } ( = 5 {{c'X u t>{))-Y x/c -,y (2.6) 

As a consequence, the semigroup densities of a semi-stable process with index 7 have the following 
property: 

Pt ( X ,y) = t-^ Pl (^lL). (2.7) 

Remark that the expression for the semigroup densities in Theorem 12.31 satisfies this property for 
7 = a/2. This remark yields the following corollary: 

Corollary 2.5. A Markov process that enjoys the time-inversion property of degree a is a semi- 
stable process of index a/2, or is in h-transform relationship with it. The converse is not true. 

Remark 2.6. Throughout the paper, we assume the semigroup densities p t (x,y) to be positive 
over the domain and regular enough to be at least twice differentiable in the space and time 
variables. This assumption is more a technical requirement for the proof than a necessity for the 
characterization of Markov processes enjoying the time-inversion property. 



3 Proof of the Theorem 
3.1 Sufficiency 

If pt (x, y) satisfies the condition l|2.4p . then from formula l|2.2p the semigroup densities q^(a,b) 
can be written as 



where 



2) H^ tf ) J « ( ^ 



s a \ t-s 



R?l = exp\p(xt^)+p[ ] + P \ -- 



at 



t a J ' \t a \ t- s 



a ( st \ 2 \ / ot \ (a 
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Using the properties of $,9,p described in Theorem 12. 3^ we obtain: 



R 



(1) = $(x,6) 

r[ 2 1 = (t-s)-^e(b) 



R 



(3) 



exp < (t — s)p(x) + p 



(t-s)- 



(t-s) 



Hence there is no separate dependence on s and t, but only on the difference t — s, which allows to 
conclude that 

and proves homogeneity for the process (t a Xi,t > 0). If in addition $(x,y) = <fr(y,x), then 



$(x,a) (t-s)^ \ (t - s)% ' (t - s) 



b b 



(t-sy. 



exp< p 



(t-sy 



(t-sy 



in which we recognize i|2.5p . 



3.2 Necessity 

For simplicity of notation, we prove the necessity of condition l|2.4p in the case a — 1. The extension 
to a > is immediate by the change of variables x i— > . Recall first the following definition 

of homogeneous functions: 



Definition 3.1. A function / : 



that satisfies for x = (x\, . . . , x n ), VA € R+, 
f(Xx) = \f*f{x) 



is called homogeneous of degree (3 £ E. If / G C 1 (M™), Euler's Homogeneous Function Theorem 
gives a necessary and sufficient condition for the function f(x) to be homogeneous: 

n A 

^, •,—/;,•) if:.r:. 

2—1 * 

Consider the function of 2n + 1 variables: Z(x, y, t) = ln(pt(x, y)), x, y S M n . From l|2.2p . Z must 
satisfy for s = t — h: 



d_ 

dt 



' 1 ,/ 6 1\ ,/& a 

t V t t \t t-h t-h t 



1 -X\-iU " 



1 



t-h t-h 



= 0. 



(3.1) 



3.2.1 The kernel $ (x, y) 

Taking derivatives with respect to bi and dj for some < i,j < n yields: 

1 1' 



d d d ^ ( b a 



dbi daj dt \t' t — h' t — h t 



= 0. 



If we set cj) (x, y, t) = -^r-^j-l (x, y, t), the latter becomes 

i ( , b X i 

7 ' "if 



t 2 (t -h)y t 



t{t - hf 



t-h 



. v. 
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1 / 1 



d(j) = 0, 



t{t-h) \t 2 {t-hf 

with the notation Vi = (g§7, ■ • ■ , £§-) T i ^2 = (gfrj • • • j af~) T ' ^ the derivative with respect to 
the time variable and = ( |, — |V For clarity, we change variables to 



b a 1 1 

z l — T ; z 2 — 7 T ; ^1 — T , ^2 



Then = 0(zi, Z2>£i + £2) and 



ti ( + 21 • Vi0 + ti 90 ) = t 2 ( + z 2 • V 2 + i 2 90 ) , 

or equivalently 

fx ( + zi ■ Vi4> + (h + t 2 ) 90 I = i 2 ( + z 2 ■ V 2 + (*i + t 2 ) 



We change variables once more: u = v = t\ + i 2 to get: 

+ 2 X -Vi0 + i> 90 = u(0 + z 2 ■ V 2 +vd(j)). (3.2) 
This gives the following proposition. 

Proposition 3.2. Equation $3.2\) is satisfied if and only if 

(f)(\x,[iy,\fit) = cj)(x,y,t). (3.3) 



Proof. As the LHS of equation (|3.2j) is independent of u, one can readily take the equivalent con- 
dition: 

+ ■ Vi0 +u 90 = 0, (3.4) 
+ z 2 - V 2 + « 90 = 0. (3.5) 

Let <?(A) = 0(Azi, z 2 ,\v) and /i(/i) = 4>(zi,fiz 2 , fxv). Equation l|3.4p implies g(X) + Xg'(X) = 0, which 
solves to g(X) = \~ 1 g(l) and hence 

0(Azi,z 2 , Xv) = j<j>(zi,z 2 ,v). 

Equation i|3.5p implies h(fi) + fih'([i) — 0, which solves to h{[i) — and hence 

0(zi,/X2 2 ,//u) = -<j){z 1 ,z 2 ,v). 



Combining the latter two equations yields i|3.3p . 

Conversely, if g is homogeneous of degree —1 and h is homogeneous of degree —1, we get: 

g(\)+\g'(\) = h(ii) + f iti(n)=0, 

which is equivalent to l|3.4p and l|3.5p and concludes the proof. □ 

The scaling property (|3.3j) implies moreover the equivalent formulation 

^'^Mttti)' (3 ' 6) 
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where <fri(x, y) = 4>(x, y, 1) must satisfy 

4>x{\x,y) = <t>i{x,\y). (3.7) 

x y — 

Under the change of variables x = y = the kernel l(x,y,t) = l(x,y,t) satisfies 

\/t \ft 

Note that the previous results remain valid for all second derivatives of l{x,y,t) with respect to Xi 
and yj, Mi,j € {1, ... , n}. By integration over Xi and yj, one can thus already make an assumption 
on the general shape of the kernel l(x,y,t): 

l(x,y,t) = k(x,y) + <pi(x,t) + ip 2 (y,t), 

where k : K 2 " — * R does not depend explicitly on time and if i, if 2 ■ R™ +1 — > M.. An equivalent 
representation for l(x,y,t) gives 

l(x,y,t) = k (jj=>-j%) +<Pi{x,t) + <p 2 (y,t), (3-8) 

with (f\{x,t) = (pi(x,t) and ip 2 (y,t) = <p 2 (y,t). 

The scaling property l|3.7p of <f>i (x, y) now translates to k(x,y) as follows: 

JLJL k {\x,y) = ^-^-k(x,\y), Vt,je {l,...,n}. (3.9) 
axi dy 3 dxi dy 3 

By integration, k(x,y) must satisfy the following condition, 

k(\x, y) = k{x, Xy) + Ci(x, A) + ( 2 {y, A), (3.10) 

along with the properties Ci( x j 1) + C2 1) = and 

(l(x,a\) + ( 2 (y,a\) = (i(x,a) + (i(ax, A) +(2(2/, A) + C 2 (A?y,a). 

Proposition 3.3. TTie functions £i(x, A) and £2(2/) A) are equivalently defined by 

(i(x,a) = d(x) - d(ax), 

( 2 (x,a) = ( 2 (x) - ( 2 {ax), (3.11) 

with a slight abuse of notation. 

Proof. The first property of the functions (i(x, A), (, 2 {y, A) implies 

Ci(z,i) = -C2(y,i) = CeR 

which is set to without loss of generality. The second property leads to 

(i(x,aX) = (i(x, a) + C,\{ax, A) + F(a, A), 

( 2 (y,a\) = C 2 (2/,A)+C2(A 2 /,a)- J F(a,A), (3.12) 

for some function F : M 2 — > R. Moreover, by deriving the latter expression involving C,i(x, A) 
successively by and A, we get 

^-^ Cl( ^ A) = |-^ Cl ^ aA) - (3 - 13) 



On the other hand, taking the derivatives with respect to A and a leads to 

d d 

C,[{x.a\) + aXQ 1 (x, aA) = x ■ \7 ax C,[{ax , X) + Q^~g^ F ( a >^' 

with ' denoting the derivative with respect to the last coordinate in A). Using (|3. 13[) . we obtain 
the following: 

d d 

g--grF(a,X) = Ci(x,aX) + aX £"{x,aX) - x ■ \I x Q[(x,aX), 

which shows that the second derivative -J^J^F(a, A) is a function of the product aA only. Hence 
with a slight abuse of notation, 

F{a,X)=F(aX)-F 1 (a)-F 2 (X), 

where we set -Fi(l) = F 2 (l) = without loss of generality. Setting respectively a = 1 and A = 1 
yields from (f3~T2|) . Fi(X) = F(X) and F 2 (X) = F(X), hence = 0. Equation (f3~T2| then becomes 

CiO,aA) -F(aX) = (i{x, a) — F(a) + (i(ax, A) — F(X), 
( 2 (y,aX)+F(aX) = (2(11, A) + F(X) + (2(Xy, a) + F(a). 

F(X) turns out to simply shift the functions C\(x,X) and C,2{y, A). It can thus be included as part 
of those functions, so without loss of generality F{X) = and the latter equations become 

CiO, aA) = Ci(z,a) + C,i{ax,X), 

C 2 (y,aA) - C 2 (y,A) + C 2 (Ay,a). (3.14) 

The latter functions can thus be further reduced to functions of one variable. With a slight abuse 
of notation, setting d{x) = Q (x, an d choosing respectively A = a \\ x \\ and a = yields 

the resulting equations l|3.1ip . □ 

This last proposition combined with (|3 . 10[) gives 

k(Xx,y) + Ci(Xx) - ( 2 (y) = k(x, Xy) + d(x) - Ca(Ay). 

Now setting $ (a;, y) = exp (k(x, y) + (i(x) — (2(y)) completes the proof of the first condition of the 
theorem, i.e. 

${Xx,y) = $(x,Xy). 

3.2.2 The function p(x) 

Going back to l|3.ip . we replace the explicit form of l(x,y,t) and use the scaling property of 
(i(x, A), C2(y, A) and k{x,y) to obtain 



d_ 

Of 



,ln i + Ci(x, Vt) - Ci(*, Vt-h)- C2 (b, -^J + C2 (a, J_ j 



ib t-h\ 1 Oi r~t~ 



V t-hj \t t(t - h) 

b 1\ (a 1 \ /a 



u^'^v ^v*-^-m 



= 0. (3.15) 



Note that the terms composed of the function k(x, y) do not depend on t and thus cancel out with 
the time derivative. 
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Recall that the variables a, b, x are independent from each other. So taking the derivative of 
p.lpp with respect to Xi gives 







*i ^^i^^i)- — — ^ [x,— ] =t 2 2 — p' 1 (x,t 2 ) 



2 ax 



i 



d 



~2 d rl 



dxi 



2 



with = ~, t% = and ' still denotes the derivative with respect to the last coordinate. Since 
the latter is valid for all ti,ts > 0, each side of the equation must be independent of t. Hence, for 
some differentiable function ipn : 1" — > K, 



d d 



1 d 



d d 



ir^Vi{^,t) = --— Vll {x) - — — Ci z,— 



dt dx 



t 2 dx 



dt dxi 



which integrates to 



<pi(x,t) = -fn(x) - Ci \x,-j=) +h 1 (x) + n(t) 



There are so far no further conditions to add on the function <pi(x, t). 

With the new result for the shape of <pi(x,t) and using the scaling property of (i(x,\) and 
£2(2/, A), we derive l|3.15p with respect to bi to obtain 



d_d_ 
dhdi 



pn+l 



Ci [b, - - C2 6, 



& 1 



t(t-h) 



= 0. 



be such that 



For convenience, let (^22 

922OM) = C2 (ri^t) - Ci (~,*) + jfn(x) - hi(x) + ip 2 2(x,t). 

The former partial derivative equation then becomes 

't 2 



¥22 



b 1 



which develops to 



a 5 rt* (b 

dbidt [h^ 11 [l 



(3.16) 



(3.17) 



with the notation z = |, ii = |. Since the LHS is independent of /i, both sides of the equation must 
cancel out. This is the case if di(p22 [z, t) is homogeneous of degree —1 and dupw (z) is homogeneous 
of degree 1, which means by integration over Zi, 



Lpn(Xz) — X 2 ipn(z) + 93 In A, <p S 



(3.18) 



Let p(z) — ipn(z), then, conditioned on showing ip — 0, we recover the second condition of the 
theorem, that is 

p(\z) = \ 2 p(z). 

3.2.3 The function 6{y) 

In order to further investigate the properties of <f22{z,t), we derive (|3. 15|) by m: 



d_d_ 

da; dt 



h 



t-h 1 t(t- h) 



f22 



+ Ci [a, t - C2 a, r 



V22 
= 



t-h t-h 



8 



using the scaling properties of (i(x,A), £ 2 (y, A) and ipn(x). Since di(f22{z,t) is homogeneous of 
degree —1, this leads to 



d_ 

dt 



t 



di^22 a, - + did a,- - dib a, - 



t 



t 



0. 



Hence, by integration over time, we obtain 



h 



t 



for some differentiable function ft. 2 



di<p22 «,- + 9iCi a, t - 9iC2 a, r = d l h 2 {a 



I, Integration with respect to <Zj yields 



1 



1 



ip 2 2 (a,t) = h 2 (a) + T 2 (t) - Ci [a, - ) + Ca [a, - 



for some real function T2(t). Furthermore the homogeneity condition on di(fi22{z,t) translates to 
dih,2{z) as follows: 

Xd l h 2 (\a) = dih 2 (a) + di( 2 (a, A) - d l C 1 (a, A), 
by the scaling properties of diQ\(x, A) and di&iy, A). Integration over then yields 

h 2 (\a) = h 2 {a) + H{\) + C 2 (a, A) - d(a, A), 

for H : M M with the properties = and F(AiA 2 ) = #(Ai) + ff(A 2 ). The function ff(A) 

is thus of the form H(\) = /31n A for some (3 £ K. Using Proposition 13.31 the latter becomes 

/*2 (Aa) + C2 ( Ao) - Cl (Aa) = /? In A + h 2 (a) + C 2 (a) - Ci (a) • 

Let 0{z) = exp (ft. 2 (z) + C 2 ( a ) — Ci( a ))- We then recover condition 3 of the theorem, i.e. 

6{\z) = \ 13 6{z). 

3.2.4 The function r(t) 

The results for ip\{x,t) and (/? 2 (y,i) summarize so far to: 



<Pi(x,t) 

<P2(V,t) 



P( x ) 
t 

p(y) 
t 



+ C 2 (y)-C 2 



+ h 1 {x)+T 1 (t) 



^ -hi(y) + T 2 (t) + h2(y). 



This implies that l(x,y,t) = \xi(p t (x,y)) has the form 
l(x,y,t) = ln$ 



( x_ y_ 



[n m + pix)+ / iy) +T 1 (t) + T2(t), 



T)(x) 



(3.19) 



where rj(x) = exp {Ci{x) — hi(x)) defines an /i-transform and can thus be neglected. It remains to 
explicitly formulate the form of the function r : K + — > R that we define by 



r(i) =r 1 {t)+T2{t) + -{n + p)\nt 



to satisfy the resulting partial differential equation <)3.15|) : 



d_ 

dt 



<p 



t(t - h) ( t 



> I - 



■ In - 

t - h t 

h 

t(t - h) 



1 t- h 



In- 



t-h 



t - h 



= 0. 



(3.20) 
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The former equation becomes after derivation, 



-fix (2 Lata + 1) + *2(21nii + 1)) = tiVfa) + {tj - ^)r'(t 2 - *i) - *2 T '(*2), 

for <i = j, < 2 = which hints that (p = 0. Assuming <p to be non-zero and setting g(t) = 

i 2 1 

— — r'(i) + 77, we get indeed 

zip 2 

h lnt 2 + *a ln*i = *i (fl(*a - *i) - ff(*i) + flf(*a)) + *a ($(*a - *i) + - ff(*a)) , 

which has no solution. So ^ = 0. Now dividing the former equation with tp = by t± and taking 
the limit ti — * leads to the differential equation 

r"(t)=0, 

which solves for r(t) = i/t, ^ € R. However, e 1 '* can always be included in an ^-transform, so i/ is 
set to and r(t) becomes trivial. 

3.2.5 The general case a > 

Combining the different factors, we obtain for the semigroup densities, 

p( x ) , p(y) 



%) exp 



f 



which simplifies to, using the scaling property of 9{y) and p(a;), 



For a > 0, we apply the change of variables cc i— > IMI^jjfir- The semigroup densities become 

Pt( x >y) =Pt (x",y°) J{y) 

where J(y) is the Jacobian of the inverse transformation, that is J(y) — a n y{" ■ ■ ■ yfi , and 
we use the slight abuse of notation x ° = (x" , . . . , x% ^ . The semigroup densities can be recast into 

where $(a;,y) = <f>^2;a,?/ = ^ satisfies condition 1, = |y«J J(y) satisfies condition 2 for 

/3 = ^7^ — n and /5(x) = p n/<*^ satisfies condition 3 of equation l|2.4p . 

4 Application to diffusions on E + 

The case of the diffusions on R+ was entirely characterized by Watanabe in [27] . It was shown that 
only Bessel processes in the wide sense (which we recall the definition below) enjoy the time-inversion 
property of degree 1. 

Definition 4.1. For some v > — 1 and c > 0, the diffusion process generated by 

Id 2 (2v + l h' c {x)\ d 



2dx 2 ' y 2x ' h c (x) J dx ^ 
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is called Bessel process in the wide sense. The function h c (x) is given by 

h c (x) = ¥Y(y + 1) (V2c x)~ v I u (V2c x), (4.2) 
where l v is the modified Bessel function. 

Remark 4.2. The Bessel process in the wide sense is in /i-transform relationship, for h = h c , with 
the Bessel process. 

We show that the result in [27] is a consequence of Theorem 12.31 which has the following one- 
dimensional formulation: 

Theorem 4.3. The Markov process (t a Xi,t > 0) on K.+ is homogeneous if and only if the semi- 
group densities of (Xt,t > 0) are of the form: 

Pt (x,y)=t-^t(^) / expj-^^ + ^j (4.3) 

for k > 0, or if it is in h-transform relationship with it. Moreover, the semigroup densities are 
related as follows: 

q ( t X) (a,b) = exp (-t £ x^ p t (a,b). (4.4) 

Proof. Theorem 12 .31 formulated on R+ gives $(x, y) = 4>{xy), p(x) = — 4j- x& , k > 0, and 0(y) = y 13 
for conditions 1, 2 and 3 to be satisfied. Equation l|4.4p is then a consequence of (|2.2p . □ 

We identify further the class of diffusion processes and provide a different proof for the result in 

Proposition 4.4. If (X t ,t > 0) is a diffusion process and (tXi,t > 0) is homogeneous and 
conservative, then both are necessarily (possibly time-scaled) Bessel processes in the wide sense. 

Proof. If (X t ,t > 0) is a diffusion process, then its infinitesimal generator has the following general 
structure: 

d d 



m{dy) ds(y) 

where m(dy) is called the speed measure and s(y) the scale measure density (see [IB])- From the 
assumptions of the proposition, (X t ,t > 0) enjoys the time-inversion property of degree 1. As a 
consequence of (|4.3p in Theorem 14.31 the speed measure is thus absolutely continuous with respect 
to the Lebesgue measure. Let u(x) — , \ , . and six) = „7^if\s - The infinitesimal generator 
can then be expressed as 

d 2 d 

where it remains to identify the functions s(y) and fJ,(y). For a fixed x > 0, let ^ be the infinitesimal 
generator of (tXi.t > 0). From equation (|4.4p . W has the following relationship with : 



which develops to 



(x) .f(b) = s(b)f"(b) + L(b) + 2x s(b) f'(b) + U(x, b)f(b). 
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For the process to be conservative, we require U (x, b) = 0, which implies no killing in the interior 
of the domain, that is 

a c//'(xb) <p(xb) k 2 2 

s (b) x , , , . + nib) x , . , . — x =0. 

^ 7 <j>{xb) PV ; </>(x6) 2 

We change variables to z = xb to obtain 

s{z/x) <j>'\z) + ^± <j>'{z) - % cf>(z) = 0. 

IE Z 

Since the latter must be valid for all x > 0, we are led to set: s(b) = ^- for a > and = ^- 
for v > — 1. This yields the following equation 

^) + ^^)-^^) = o. 

The general solution (non-singular at and up to a constant factor) is expressed through the 
modified Bessel function of the first kind as follows: 

Gathering the different factors in l|4.3p leads to 

(?)"" M-t(t + ? 

where TV is a normalization factor. The additional condition that lim Pt(x, y) — S(x — y) implies 

j3 = 2v + 1 and k = — , which leads to the semigroup densities of a time-scaled Bessel process of 
dimension ^: 



x 2 + y 2 
2<jH 



pt{x ' v)= ik (IT ^S) exp 

The infinitesimal generator for (tXi , t > 0) is given by 



2 <% 2 [26 J 56' 

where one recognizes expression l|4.ip for h c (b) = 4>(xb) and a time-scale t >— > cr 2 i. □ 

Proposition |4]4] smoothly extends to any power a > of the Bessel process through the mapping 
x i— > ||£|| a TT§T|. In particular, it is worth remarking that the case a = 2 gives rise to squares of 
Bessel processes, which leads to the following result: 

Proposition 4.5. // {X t ,t > 0) is a diffusion process and {t 2 Xi,t > 0) is homogeneous and 
conservative, then both are necessarily (possibly time-scaled) squares of Bessel processes in the wide 
sense. 

5 Examples 

5.1 Generalized Dunkl processes and Jacobi-Dunkl processes 

5.1.1 Multidimensional Dunkl processes 

We briefly review the construction of the Dunkl process in 1™ (see [211 26j). 
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Definition 5.1. The Dunkl process in M™ is the Markov cadlag process with infinitesimal generator 

i=l 

where T i: 1 < i < n, is a one-dimensional differential-difference operator defined for u g C 1 (R n ) by 

(•,•) is the usual scalar product, i? is a root system in W L and i? + a positive subsystem. k(a) is 
a non-negative multiplicity function defined on R and invariant by the finite reflection group W 
associated with R. a a is the reflection operator with respect to the hyperplane H a orthogonal to 
a such that a a x = x — (a 1 x}a and for convenience {a, a) — 2 (see [HUH])- 

A result obtained by M. Rosier [25J yields the semigroup densities as follows: 

M fn . 1 _ ( \*\ 2 + \v?\ n. ( x v 



« = D > U' v! » Wfc(2/) (5 - 3) 



where Dk{x,y) > is the Dunkl kernel, uJk{y) = j~J | (a, y) | 2fc ( Q ) the weight function which is 

k{a) and Ck = e ~ uik(x)dx. 

aeR+ 

Following a thorough study of the properties of the one-dimensional Dunkl process in [12], it 
was remarked in [13] that the Dunkl process in K" enjoys the time-inversion property of degree 1. 
Considering that the Dunkl kernel satisfies 

D k (x,y) = Dk{y,x) and D k (fix,y) = D fc (x, fiy), (5.4) 

the proof is straightforward with 

$(x,y) = D k (x,y), 6{y)=u k {y), p{x) = (5.5) 

By equation l|2.5p . the semigroup densities of the time- inverted process is even in /i-transform 
relationship with the semigroup densities of the original Dunkl process: 

(*, (a ,5 ) = ^4ex P (- M! t ) pW(a,6). (5.6) 
D k (x,a) \ 2 J 

5.1.2 Generalized Dunkl processes 

In an attempt to generalize the Dunkl process, we extend the definition of the infinitesimal generator 

to 

f{x) = 2 A/(S) + E fc ^ (g , a) + L ( ) (:E>a) 2 (5-7) 

where A is the usual Laplacian, / € C 2 (R") and X(a) is a non-negative multiplicity function defined 
on R and invariant by the finite reflection group W, similarly as k(a). We retrieve the Dunkl process 
for X(a) = k(a). Note that these processes are no longer martingales for X(a) ^ k(a). 

The one-dimensional case was introduced in [12], where the semigroup densities were explicitly 
derived, 

P^\x,y) = /-exp (-^f) D KX ) (5.8) 
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with the generalized Dunkl kernel, [y = k 



V^ 2 +4A), 



(5.9) 



for z 



. From (|5,4p . the generalized Dunkl kernel satisfies 



D k ,x(x, y) = D k ,\(y, x) and D k ,\{/J,x, y) = D kt \{x, ny), 



(5.10) 



which readily implies that the generalized Dunkl process also enjoys the time-inversion property of 
degree 1. The semigroup densities were derived as an application of the skew-product representation 
of the generalized Dunkl process (X t ,t > 0) in terms of its absolute value (a Bessel process) and 
an independent Poisson process TV, 



W. 



x t ^ \x t \(-iy 



(5.11) 



where A t = j* ||. 

In the n-dimensional case, the application of the skew-product representation derived by Chy- 
biryakov |6j shows that the generalized Dunkl process enjoys time-inversion for some specific root 
systems. We first recall one of the main results of [6]. 

Proposition 5.2. Let (X t ,t > 0) be the generalized Dunkl process generated by Q5. 7\) . with (X^ ,t > 
0) its radial part, i.e. the process confined to a Weyl chamber. Let R+ = {«i, . . . , a{\ for some I € N 
be the corresponding positive root system and let (NJ:,t > 0), i = 1, ... ,1 be independent Poisson 
processes of respective intensities A(«i). Then X t may be represented as Y t l , which is defined by 
induction as follows: 



Y t " = Xf and Yl = a a 



y: 



1, 



,1, 



where A\ 



ds 



o (Y 



The proof follows the argument in [gj, while replacing k(eti) by \(ai) appropriately. 

From now on, let R+ = {a±, . . . , a{\ for some I < n be an orthogonal positive root system, that 
is (ai, ctj) = 2Sij. For this particular root system, one can show that the generalized Dunkl process 
enjoys time- inversion of degree 1. We first prove the following absolute continuity relation: 

Lemma 5.3. Let (X™ ,t > 0) be the radial Dunkl process with infinitesimal generator 



(5.12) 



Fixed v e {1,...,/}. Let k'{a) be another coefficient function on the root system i?+ such that 
k'(a.v) > k(a u ) and k'(ai) — k(cti) for i ^ v. Then, denoting Pi the law of the radial Dunkl 
process X^ starting from x, we have 



p(fc') 



k'{a v )-k{a v ) 



exp 



a v , x) 



(k(a v ) 



ds 



(a v ,X, 



W\2 



(5.13) 



. p(fe)| 



Proof. Let ko(a) be a coefficient such that ko{a v ) = | for some fixed i S {1, . 
following martingale decomposition (see |14j): 



.,/}. X™ has the 



X? = x 



i 

j>(«i) 

i=l 



ds 
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where B^ is a (P^ 01 , JF^-Brownian motion. Consider the local martingale 



L t (fe,) = exp( (fc'K)-~ 



(a v ,dB, 



(feo) 



i) 2 







for some coefficient function k'(a) such that k'(a v ) > | and fc'(aj) = fco(c*i) for i 7^ ;a The Ito 
formula for In ((a v , X^)) combined with the orthogonality of the roots yields 



L 



(k>) 



[otv, X. 



(a u ,x) 



exp 



(fc'K)-i) 2 r* 



ds 



{a v ,X f 



W\2 
>-"-s / 



Define the new law By the Girsanov theorem, 



B 



(fc') 



Si 



is a (Pi fc ; , JF t )-Brownian motion and hence, 



Xf = x + Bp 



is a radial Dunkl process under (P?' 3 ,^)- 

Define fc as another coefficient on the root system that satisfies the conditions enunciated in the 
lemma. The absolute continuity relation is then a consequence of the successive application of the 
latter result to the indices k' and fc. □ 



k'{a v ) - - 



Jo 



ds 



{<*u,xy) 



ds 



Now as an application of Proposition I5.2( we prove the following: 

Proposition 5.4. Let R+ = {a\, . . . ,a{\ be an orthogonal positive root system for I < n. Then 
the generalized Dunkl process {Xt,t > 0) generated by \5. 7|) enjoys the time-inversion property of 
degree 1. 

Proof. Using orthogonality of the roots, remark that 

(ai,o-jx) 2 = (oti,x - (a j ,x)a j ) 2 = (a,, a;) 2 , 

which implies in particular 

so that the inductive representation of X t in Proposition 15.21 becomes 



*t=n o x, 



for A\ = 



(xY,<*i) 



The radial part of a Dunkl process enjoys the time- inversion property of degree 1. We need to show 
that the semigroup densities of the generalized Dunkl process are related to the semigroup densities 
of its radial parts. For / 6 C 2 (R n ), 

E x[/0?)] = E K [/(F t l_1 )l {Ar ^ is cvcn} ] + E x [f(a ai Y^ )l^ N i^ is odd }]. 
Since P(iVJ is even) = | (1 + exp(— 2A(ojj)u)), we obtain 



E x [/(*?)] = 

+ E x 



/(F/- 1 ) -(l + exp(-2A(a,)^)) 



fi^Yt 1 ) -(l-exp(-2A(a i K)) 
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The expectation [/(X t )] can thus be evaluated by induction on i £ {1, . . . , I}. It follows that the 
semigroup densities of X t can be expressed as the product of the semigroup densities of its radial 
parts times a function involving expectations of the form 



cxp (-2A(a„)A£) 



xf = y 



for v G {1, . . . , I}. From Lemma OO 



E (fc) 



exp (-2\(a v )A v t ) 



xr = y 



p[ k \x,y) S ' {a v ,y)\ k{ai,) k{oi "> 



p { t k) (x,y) \(a v ,x) 



where fc'(o;,) = k(a>i) for i ^ v and k'(a u ) = | + J (k(a,y) — i) +4X(a^). The form of the 



semigroup densities in (|5,3p implies that the expectation is a ratio of Dunkl kernels, 

n,, ( x y\ tm, ( v \ / 1„ v \ \ ( k - k ')( a ») 
Uk ' [Vt'Vi) Wk ' \~t) / < a "' ft' \ 



E« 



-2\{a v )A- t 



xf = y 



Ck 



which reduces to 



E (fe) e -2A(a„)A t " 



Xf = y 



c k Dk ' {ft 1 7i) 



Ck' D , 



JL\i JL' 
av > y /i ){a '" ft' 



(fe'-fc)(av) 



by definition of k'(a). The conditional expectation thus satisfies condition 1 of Theorem 
consequence, the conditions of Theorem 12.31 are satisfied for 



As 



®(x,y) = D kt \{x,y), 6(y)=ujk(y), p(x) 



\ x \ 



where Dk^xix^) is a generalized Dunkl kernel given explicitly in terms of radial Dunkl kernels and 
satisfying equivalent conditions (see (|5.4p ). □ 



5.1.3 Jacobi-Dunkl processes 

Gallardo et al. j5j derived the Jacobi-Dunkl process as the hyperbolic analog of the one-dimensional 
Dunkl process. It is defined as the process generated by 



J[ ' dx 2 A{x) dx dx\A{x))\ 2 ) 



(5.14) 



where A(x) = (sinh 2 (x)) 2 (cosh 2 (x)) 2 . From the expression of the semigroup densities de- 
veloped in |5|, this process does not enjoy the time-inversion property. Its radial part however 
corresponds to the Jacobi process of index (a,/3) on R+ (see [l6t I17|). The infinitesimal generator 
of the Jacobi process, expressed by 



is in /i-transform relationship with the Laplacian operator for h(x) — y/A(x). Since the Brownian 
motion enjoys the time- inversion property of degree 1, so does the Jacobi process by Theorem [ 
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5.2 Matrix-valued processes 



5.2.1 Eigenvalue processes 

Dyson in [11] described the eigenvalues of a Hermitian Brownian motion as the joint evolution of 
independent Brownian motions conditioned never to collide (see also [15] and [4]). It was further 
remarked that the process version of the Gaussian orthogonal ensemble does not admit such a 
representation for its eigenvalues. This work was extended by Konig and O'Connell in [19] to the 
process version of the Laguerre ensemble, denominated the Laguerre process and defined as follows: 

Definition 5.5. Let Bt be an n x m matrix with independent standard complex Brownian entries. 
The Laguerre process is the matrix-valued process defined by {X t = B' t B t ,t > 0}, where B' t is the 
transpose of B t . 

From [T5J, the eigenvalues of the Laguerre process evolve like m independent squared Bessel pro- 
cesses conditioned never to collide. No such representation however exists for the case where the 
entries of B t are real Brownian motions, i.e. the Wishart process considered by Bru in [3]. 

The main result of [19] is that both of the above mentioned eigenvalues processes (in the complex 
Brownian case) can be obtained as the /i-transform of processes with m independent components. 
The joint eigenvalues process is thus in /i-transform relationship with a process that enjoys the time 
inversion property of degree 1 in the case of the m-dimensional Brownian motion and degree 2 in 
the case of the m-dimensional squared Bessel process, as made explicit in the following proposition: 

Proposition 5.6. Letpt(xi, yi) (i = 1, . . . , m) be the semigroup densities of squared Bessel processes 
(respectively Brownian motions), and let 

m 

h{x) = Y[(Xj - Xi ) (5.16) 

i<j 

for x — (#1, . . . , x rn ). Then, the semigroup densities of the joint eigenvalues process of the Laguerre 
process (respectively the Hermitian Brownian motion) are given by 

Pt(x,y) = J^j det (Pt( x iiyj))™j=i ( 5 - 17 ) 

m 

with respect to the Lebesgue measure dy = J^J dy.j . 

3 = 1 

It follows immediately by Theorem 12.31 that the eigenvalues processes enjoy the time-inversion 
property. Moreover by Proposition 12.21 they yield the same process under time-inversion as the 
m-dimensional Brownian motion or the m-dimensional squared Bessel process respectively. 



5.2.2 Wishart processes 

The Wishart process WIS(<5, tl m , jx), introduced by Bru in [T], is a continuous Markov process 
taking values in the space of real symmetric positive definite m x m matrices It is solution to 
the following stochastic differential equation: 

dX t = y/X~ t dB t + dB' t y/x~ t + SI m dt, X = x, (5.18) 

where Bt is an m x m matrix with Brownian entries and I m the identity matrix. Further results 
have been obtained in [2] and [3]. In [7], among other major findings about the Wishart process, 
the transition probability densities expressed with respect to the Lebesgue measure dy = Y\_(^Vij) 

were derived in terms of generalized Bessel functions (we refer to the appendix for the definition) : 



(5.19) 
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for x,y S S'jjj and S > m — 1. From the shape of its densities, the Wishart process was stated 
in [13] as an example of Markov processes enjoying the time-inversion property of degree 2. The 
hypothesis of Theorem 12.31 is indeed satisfied for n = ^m(m + 1) and 



®{x,y) = (det(x) det(y))~ 



6(y)=l-(det(y))^ 



-1 - 

~ li 



1, 



> p( x ) = —2 Tr ( x )- 



(5.20) 



Next we use a skew-product representation, as for the Dunkl process, to elaborate on the Wishart 
process and derive a matrix-valued process with jumps. The skew-product representation allows 
the expression of the semigroup densities in terms of the Wishart transition probability densities. 

Definition 5.7. Let (N^\t > 0) be a Poisson process with intensity A. Let {X t ,t > 0) be a 
Wishart process W\&(8,tl m ,\x) independent of the Poisson process. The skew- Wishart process 

(Xf , t > 0) is defined through the skew-product 



X 



x t {-if** 



(A) 



(5.21) 



where A t 



Tr(X- l )da. 



Proposition 5.8. The transition probability densities of the skew-Wishart process are related to 
the semigroup densities p t (x,y) of the Wishart process X t as follows 



Pt(x,\y\) I ^ {yeS +}\ \ l+ \T 



-xy 



At 2 



(5.22) 



for v = i=™=l, t/ = V^+4A and \y\ = y(l {yeS + } - l {s6S - } ). 

Proof. Let (Pt)t>o be the semigroup of the skew-Wishart process. For x > and / € C c (M m (M)), 
P*/(z) = E, 



/(X t (A) ) 



E r 



f( X t) 1 {N M is evcn} 



f\~Xt) lrurM • 

J \ <•! \N\ ' IS o 



dd} 



With P(Ni X) is even) = |(1 + exp(-2Au)), we have 



Ptf(x) 



E, 



/(X t ) ~(l + exp(-2AA t )) 



/(-X t ) ~(l-exp(-2AA t )) 



(5.23) 



(5' —7TL— 1 



Let Q^' 3 with i/ 

with = $=sf=I 
the probability laws are related as follows: 



denote the probability law of a Wishart process WIS(S',tI, \x), and 
the probability law of X t . According to Theorem 1.2 (Remark 2.3) in [7], 



Qx \j? t — 
from which we deduce 



/dot X t 
\ detx 



exp 



Tr(X-')ds ) 



>') 



(x,y) ( dety 



p["\x,y) Vdet: 



Q 



(") 



exp 



TriX-^ds 



Xi 
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Thus, from the expression of the semigroup densities in (|5 . 19[) . we have 



exp (~ 2A ^ ^(X'^ds 



x t = y 



Combining the latter with l|5.23p yields the semigroup densities for the skew-Wishart process. □ 

The skew-Wishart is an example of matrix-valued process with jumps that enjoys the time- 
inversion property of degree 2. Indeed, by setting 

*(x,y) = (det^detOi/l))"* |l {ygs + } i(t+t')(^) 
+ 1 {yes-}2 ~ 



-xy 



0( y ) = ±-(det(\y\)y , p(x) = ~TrQx\), 



the conditions of Theorem 12.31 are satisfied for a = 2 



(5.24) 



A Generalized hypergeometric functions 

Using the notation in Muirhead [5T| , hypergeometric functions of matrix arguments are defined for 
a real symmetric m x m matrix X , a, € C and 6j S C\{0, i, 1, ... , 222 j^} by 



,Fq(ai, . . ■ , a p ; 6i, . ■ ■ , b q ; X) = }^ 



k=0 K 



(ai) K ■ ■ ■ (a p ) K C K {X) 



(A.l) 



where the second summation is over all partitions k = (ki, . . . , k m ), ki > ■■■ > k m > 0, of 
k = YliLi k\ = k\\ - ■ ■ k m \ and the generalized Pochhammer symbols are given by 



= n 



i - 1 



(a)k = a(a + 1) •• • (a + k — 1), (a) = 1. 



C K (X) is the zonal polynomial corresponding to k, which is a symmetric, homogeneous polynomial 
of degree k in the eigenvalues of X that satisfies 



c k {yx) = c k (VyxVy) 



(A.2) 



for some Y € The function p F q (ai, . . . , a p ; b\, . . . , b q ; YX) thus makes sense. Finally, we define 
the generalized modified Bessel function by 



r,„(- + ^i; 

where the generalized gamma function is given as a product of the usual gamma functions, 

i - r 



T m (a) = 7r" 



(A.3) 



(A.4) 



for Re(a) > ^j^- Note that the generalized modified Bessel for m = 1 relates to the usual modified 
Bessel function J„(x) by I u (x) = I v {2^/x) (see [23J). 
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